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ABSTRACT 


A searchlight type sonar is one of the systems that 
small navies use to counteract the danger which submarines 
present to their lines of supply and transport. 

In this paper, a standard search pattern for this 
type of sonar is compared with search patterns which are 
based on a consideration of the tactical value of detecting 
a submarine as a function of the relative location of the 
Submarine. 

The results of the comparison suggest that is possible 
to increase the effectiveness of a Secanech ight. type sonar 
by uSing a search pattern in which the Sweep time allocated 


to a search sector is based on the sectors tactical importance. 
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ieee Nt RODUCTION 


Antisubmarine Warfare (ASW) training constitutes a 
major part of many navies peacetime training. One reason 
for this is that the ability to transport troops and supplies 
by sea during wartime could be a decisive factor in determining 
the outcome of the conflict. 

Several ways exist for navies to improve their ASW 
capabilities other than to obtain new ASW sensors and weapons 
systems; one of these is to develop improved tactics for use 
with their existing systems. 

An important ASW sensor system is shipboard mounted active 
sonar. Two types of active shipboard sonar systems which are 
in use today in some small navies are the searchlight type 
sonar system and the scanning type sonar system. In the 
former, a sound pulse is transmitted into and echoes are 
received from a narrow sector in a given direction. The 
area surrounding the ship must be searched in steps, sector 
by sector. In the scanning sonar, a sound pulse is transmit- 
ted and echoes are received in an omnidirectional way. The 
area surrounding the ship is searched on each pulse. 

For the searchlight type sonar, the operator has control 
of the pattern to be followed in searching the sectors 
surrounding the ship. In this paper, an investigation of 
patterns for a searchlight type sonar which tend to optimize 
detection in certain prefered sectors is made by using com- 


puter simulation techniques. 





Mie Naw oRe OF THE PROBLEM 


The problem of finding an optimum search pattern for a 
searchlight type detection system was addressed by Koopman 
(Ref. 1). As a measure of effectiveness, Koopman used W, 
the effective sweep width, and he concluded that when the 
instantaneous probability of detection during dt is jdt, 
assuming no previous detection, and i = Ay (X)/ 2, where 
Ay (x) decreases with increasing r, the search pattern which 
Maximizes the probability of detection for an infinite 
straight line encounter model consists in fixing the line of 
sight (or sonar beam) directly along the axis of abscissas, 
Gividing the time equally between the right and the left 
axis. He made the following remark "It would be misleading 
to conclude that scannings should always be confined to the 
beam. In most cases, it is imperative to detect the target 
Callin... 

In this paper, some search patterns, which consider the 
tactical value of detecting a sukmarine as a function of its 
relative location are investigated. 

A standard pattern for searching with a searchlight sonar 
1S tememert on the starboard beam and, after the first emission 
or "ping" is sent out, train the projector ten degrees toward 
the bow, send the second "ping", again train ten degrees 
toward the bow, etc., until the projector is aimed directly 
toward the bow. From that position, the projector is then 


trained to the port beam and the above procedure is repeated. 
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The time between pings is determined by the range scale 
selected by the operator which is usually the maximum range 
for the equipment. 

In order to compare this search pattern with other 
patterns, an area around the ship was divided into sectors 
as shown in Fig. 1. The semicircular boundary of the area is 
at a distance are from the ship and it is called the maximum 
range of the sonar. It will be defined more explicitly on 
page 12. 

In each sector, for the purpose of illustrating the 
approach being considered here, a somewhat arbitrary relative 
value for the detection of a submarine in the sector was assign- 
ed. The value was intended to indicate the effectiveness of a 
detection in the sector in preventing a submarine from com- 
pleting its mission. The rationale supporting the relative 
values waS in part as follows: a subinarine detected in Sec- 
tor B iS ina favorable position to be attacked, while a detec- 
tion in Sector C provides only a minimum time to attack. A 
detection in Sector A is intermediate in value to detections 
in Sector B and Sector C. 

To make the comparison between the standard pattern and 
the pattern to be conSidered here, no search effort was assign- 
ed aft of the beam, also a minimum range of detection for a 
Sonar will be assumed. These considerations are the basis for 
defimmingemseector |) in which no detection can occur. If an appli- 
cation of the approach considered here were to be made, the 
numbers of sectors, their limits and values attached to them 
should be determined by specific tactical considerations. 
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PIGUEE 1 , 


SEARCH SECTORS WITH THEIR RELATIVE VALUES 





A. MATHEMATICAL MODEL 
The values which have been assigned to the different 
sectors have the following order Va7VarVe and they will 


now be specified as being all positive. 


A negative V, was assigned to the event the target 


D 
reaches Sector D, that is, the event the target enters the 
sonar search area but is not detected. 

The problem then is to try to maximize the expected value 
of the value of detection by choice of the sweep pattern. 

To obtain the detection probability on a single ping, 


a model which is outlined by Urick (Ref. 2) was used. 


The sonar equation, as given by Urick is 


(1) OmageS/Ne = Skene 20) + TS ~ (NL = DI) 


where S/N, is the signal-to-noise ratio, TL is the trans-~ 
mission loss, TS is the target strength, SL is the source 
level, NL is the noise level and DI is the receiving directi- 
Vii Gyeenaex . 

In simple sonar detection models, S/N, is related to 
the probability of detection Pa and the false alarm proba- 
bility Pee If a false alarm probability is specified, 
then the probability of detection becomes a function of the 
Signal-to-noise ratio alone. By using this function and 
the sonar equation, the probability of detection can be 
expressed as a function of the transmission loss, the source 
level, the target strength, the noise level and the directi- 


Val Ey anGlex . 





The transmission loss is a function of the range of the 
target and for a given range r in yards it is assumed here 
to be given by 


(2) iS DO os <¢ Seem r- 


where the first term represents a spherical spreading loss 
and in the second termaisS in units of decibels per kiloyard 
and the second term represents absorption loss. 

Given such a relation between TL and range r of a target 
and values for SL, TS, NL and DI and a relation between the 
probability of detection Pa and the signal-to-noise ratio 
S/Not the probability of detection Pa Or ans ing le. pincecan 
be expressed as a function of the range to the target, that 
is the relation Pa = Pq (x) can be found. 

To determine such a relation, the characteristics assumed 
for the searchlight type sonar were the following: an output 
of 1000 watts of acoustic power, a beam width of 16 degrees, 
a pulse length of 0.1 sec at a frequency of 10 kHz. and a 
receiver bandwidth of 500 Hz. The sonar was assumed to be 
hull mounted on a destroyer traveling at 18 knots. The 
angular width of the beam was defined as the angle between 
the 38a. down rays. 


The sonar equation is often written as 


(3) TT, = 1/2 (SL + TS - NL + DI - DT) 





where DT = 10 log eae op and (S75) op is the signal-to- 
noise ratio required to achieve some operating point, that 
is, to achieve a particular pair of values of (DerPg)- The 
Significance of the equation in this form is that it relates 
the values of TL, SL, TS, NL and DI to a particular signal- 
to-noise ratio and, hence, when a relation between S/N. and 


Pe and P, have been specified, to a particular operating 


d 
pOine., For example, if an operating point (PerP,) is desired 
and a relation between S/N o+ Per and Pg and values of SL, 

TS, NL and DI are specified, then the value of TL which will 
give the desired operating point can be found. The value 

or Pa at the operating point is often said to be the value 
required in order that the targetwill be "just detectable". 

In this paper, it will be assumed that if Pa is less 
than .1, then it is in effect zero. This implies that, 
for a given Per there iS a maximum range of detection a: 
if Pg iS a non-increasing function of the range, which will 
be the case for the model used here. 

A value for Pe was chosen by specifying that the proba- 
bility of more than two false alarms in one hour was .1l. 
This implies the value of Pr is l Goie. This is shown in 
Appendix A. 


The Source Level can be expressed as (Ref. 1, p. 63) 


eee ioe + 10 log P+ Dia. 


For an acoustic power of 1000 w anda DIn = 20 db for the 
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central ray of the beam, SL = 118.5 db where the reference 


le eee Prom reference ay Fig 11.10, 


level is 0.64(10)~ 
the self noise level NL for a destroyer at 18 knots is -40 
db ateeZ5— kHz 

A noise level slope of -6 db/octave will be assumed in 
this paper. To compute the value of NL at 10 kHz, consider 


the following argument. Let x represent the number of octaves 


between 10 kHz and 25 kHz, then 
ee — 
Ty = f/f. eee 
and x = 1.32. Hence at 10 kHz 
Meee — 405—(-6) (17382) = = 32.08 db. 


The target will be assumed to be a submarine at bow 
aspect with target strength TS = + 10 db (Ref. 1, Table 9.2). 
The following model was adopted to relate signal-to-noise 


Latero: £O Pr and Pa (Ref. 4) 


(4) Dawe ead) 
(5) Pe = 1 => (re) 

Z 
(6) d = wt (S/N) 


where w is the bandwidth, t is the pulse length and N = WN os 
For a given value of Pr which in this case is he (UO) Vy 


is determined by equation (5), and with v, and given values 


ie 


Of Weand “Gy Dg became a function of S/N, through equation 


(6) Mand (4). 
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This model corresponds to that for a narrow band gaussian 
Signal in gaussian white noise where S/N<<l. The function 
Pa = Pg (r) could be defined by first computing values of 
TT, for various values for r by using equation (2) and then 
using equation (1) to compute values of 10 log S/N. With 
these values and equation (6) values of d could be computed 


and with the value of v, obtained from equation (5), the 


e 
values of Pg which corresponded to the various values of r 
could be determined by using equation (4). 

Because of the way the function was to be used, particular 
values of Py were first specified and then values of r which 
corresponded to them were determined. These values, along 
with the corresponding values of TL are listed in Table l. 

The value adopted for es: is that determined by the 
operating point 5 RON Ee sel ee 

Since the sonar equipment of the type being considered 
here is usually limited in its detection capability for 
targets at very close range, a minimum detection range 
of 200 yards was adopted in the model. 

To simplify the problem, it was assumed that the target 


waS a point target and that if a target is in the beam for 


Meoemac,, the probability of detecting it is 
m 


e a _ TI = 
P {Detection} = 1 al ee Pq (rs) } 


where both m and r, are determined by the relative track 
of the target through the area scanned by the sonar beam. 

As a further simplification, the targets track was assumed 
to be parallel to the ships track, so that a target will enter 


eZ 





2 eS OC Ge Betas wm 


0 

OF 
OF 
OF 
O: 
Or 
On 
0. 





Table 1. Values of the transmission loss and range 

for particular values of the Probability of detection 
and a false alarm probability such that the probability 
Of more than two false alarms in a one hour DOr LOG arc 


equal’ to .1. 


1S 





the semicircular area swept by the beam if its lateral range 
x is such that es | <a 

In order to determine the number of scans m on a target 
during its straight line encounter with the searcher, the 
sonar position in its search cycle at the moment the target 
crosses the semicircular boundary must be specified. In 
addition, the targets lateral range, sonar scanning period 
and the sweep pattern must be specified. It is reasonable 
to assume that given the target enters the search area, its 
lateral range is a random variable which is uniformly dis- 
tributed between sess and sae and that the target crosses 
the line which is tangent to the semicircular area swept by 
the beam and parallel with the x-axis at a time which is 
uniformally distributed over the time for a complete scan 
evele. 

Even with the above assumptions, the problem of determin- 
ing the average of the value for a pattern is difficult to 
handle analytically. However, a Monte Carlo simulation can 
provide a satisfactory solution to this problem. 

In this investigation, a computer was used to do a 
Monte Carlo simulation in order to eStimate the expected 
value of the value for some patterns. 


Tis Simulation and its reSults are discussed in the 


remainder of the paper. 
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Per eeeomocecle rT PON OF THE SIMULATION 


A search pattern which might be a realizable optimum 
one under the assumed conditions can be described as follows: 
start with one normal sweep beam to beam, and then allocate 
one partial sweep between the rays which bound Sector B. 
This will be called a 2/1 search pattern since Sector B is 
searched twice for each search of Sector A. Using this 
terminology, the standard search pattern would be called a 
1/1 search pattern. By allocating additional partial sweeps 
to Sector B, the search patterns 3/1, 4/1, etc., are generated. 

It has been assumed that a submarine can enter the sonar 
area of detection at any position of the projector and at 
any lateral range. To Simulate this, a target was generated 
at random on a line parallel to the x-axis and tangent to 
the semicircular area swept by the beam each time the projector 
was advanced to a new position. 

To reduce starting transient effects and to allow each 
target an opportunity to be detected, the search process 
was not started until the first target generated had advanced 
approximately 1750 yards and the last target generated was 
allowed to travel the same distance before the search process 
stopped. 

The simulation was run for the search patterns 1/1, 2/1. 
3/1, 4/1, and 5/1 with 500 targets generated for each pattern. 
Ten different runs of each pattern were made in order to 


obtain a statistically adequate sample. 
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The sectors chosen and the detection values assigned to 
them were as follows, where bearings are relative to the 
heading of the ship: Sector A: between bearings 040° and 
090° and between 270° and 320° and between ranges of 800 
and 1800 yards. Value: 2. Sector B: between bearings 
320° and 040° and between ranges of 800 and 1800 yards. 
Value: 10. Sector C: between bearings 270° and 090° and 
between ranges of 200 and 800 yards. Value: 1. Sector D: 
between bearings 270° and 090° and between ranges of 0 to 
200 yards and all the area aft of the x-axis. Value: -l. 

The geometry used for the simulation is shown in Fig. 2. 
The ship is at the origin of a system of rectangular coordi- 
nates which 1s moving in the positive direction of the y-axis 
at a speed of 18 knots. 

As only positive values of x need to be considered because 
of the symmetry of the sectors and sweep patterns, the area 
of interest for the simulated search is the positive quad- 
rant and targets were uniformly generated on the positive 
section of the line parallel to the x-axis. This section 
of the line is 1800 yards long and it is 1800 yards from 
the ship. The value 1800 was used since it is a little 
larger than the range corresponding to Pes 0.1 which was 
arbitrarily said to be equivalent to Pa = OF 

The targets moved in the negative direction of the y- 
axis and parallel to it at a speed of 10 knots, so that their 
relative speed during an encounter was 28 knots. Both ships 


kept their courses and speed constant. 
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The name of the variables and data used are given in 
Appendix B and the schematic flow chart is shown in Fig. 3. 
The matrix TAR and vector ITAR were used to keep the 

parameters of each target generated. The coordinates of 
the position of the target were kept in TAR. The condition 
of the target was said to be active if it was in the detec- 
tion area and not detected and passive if it was in the 
detection area and detected or if it had crossed into 
Sector D undetected. 

When the search process starts, the first ping is on 

the starboard beam, for each ping after the first ping, the 
projector is trained 10 degrees and the target advanced 
downward 35 yards which corresponds to the relative distance 
traveled by the target in a time interval equal to that 
necessary for a sound pulse to go to 1800 yards and return. 

A sound speed of 4800 ft/sec was used. 

On each ping, the angular positions of the active targets 
were computed and compared in order to determine those, if 
any, which were inside the sonar beam. A beam width of 
16 degrees was used. All of those targets which were out- 
Side of the beam could not be detected and they were advanced 
downward. If a target was inside the beam, its range to the 
ship was computed. If this range were less than the minimum 
range of 200 yards, no detection was possible and the target 
was then assigned to Sector D. 

For targets whose ranges were greater than the minimun, 
the transmission loss to the target was computed and the 


probability of detection was obtained. If this probability 
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was less than 0.1, no detection occurred and the target 
was advanced. 

To determine if a detection event occurred, a random 
number, uniformly distributed (0,1), was generated and if 
the random number was less than the probability of detection 
a detection occurred. Otherwise a detection did not occur. 
Detections were assigned to the sectors in which they 
OCeirered . 

After the above process was completed, the projector 
was trained, a new ping emitted and the process was started 
again. 

Whenever the projector reached the position in which it 
was trained toward the bow of the ship, the search on the 
port side was simulated. This was done by executing the 
generation and advance of the targets in the positive quad- 
rant for the numbers of pings required to cover the port 
side for the pattern being used. Detection of targets could 
not be made during the period necessary to generate these 
pings. The above procedure was possible because of the port/ 
starboard symmetry of the problem. 

The targets which crossed the area without being detected 
were assigned to Sector D. 

The search pattern execution was controlled by the part 
of the program labelled Control]. of the Search Sweep Pattern. 
It essentially counted the times the projector went into 
each sector and produced the desired sweep ratio. This 
section also controlled the simulated search on the port 


Side. 
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When the five search patterns had been realized, a new 
run for all the sweep patterns was executed. 

The detection values obtained for each sweep pattern 
are kept in the matrix STORE. This matrix and the average 
of the value and standard deviation of the value for each 
pattern was the output of the simulation. 

The outputs of a given search pattern for the ten 
aifferent runs can be considered to be samples from the 
population of values for that sweep pattern. 

The Mann-Whitney test was performed with the data in 


the matrix STORE and the result was printed. 
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IV. CONCLUSIONS 


The averages of the values of detection for the five 
sweep patterns investigated and for the assigned values 
suggested in some cases that the expected value of the value 
of detection is different than that of the 1/1 standard 
pattern. To test this conjecture, the Mann-Whitney test 
was used. This test was used because it is nonparametric, 
the distribution of the value of detection was not known, 
and the sample size was not large enough to assume a normal 
Giks Glade Om, 

The Mann-Whitney test (Ref. 3) with the assumption that 
if there is a difference between population distribution 
functions, then that difference is a difference in the means 
of the distributions allows one to make the following hypo- 


thesis: 


eel oe Bin 1) 


H E{X(1/1)} < Efy(n/1) } 


where X(1/1) is the random variable representing a value 
from the population of values of detection with the (1/1) 
Standard search pattern and Y(n/l) is the random variable 
representing a value from the population of values of detec- 
tion obtained with a (n/l) search pattern. 
If the null hypothesis is accepted, that is if the 

expected value of X is greater than or equal to the expected 
value of Y, this means that at the level of significance 
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chosen, there is not a significant improvement in the value 
of detection with that pattern. If the null hypothesis is 
rejected, this means that there is a significant improve- 
ment by using that sweep pattern. 

For the test,a 90% confidence level was chosen. 


The Mann-Whitney test results are given in Table II 


below. 
TABLE II 

RESULT OF MANN-WHITNEY TEST 
VALUE OF T STATISTIC 50.00 
FOR SWEEP PATTERN 1/1 ACCEPT Hy 
Wau OF T STATISTIC 43a 50 
FOR SWEEP PATTERN 2/1 ACCEPT H) 
VALUE OF T STATISTIC 36.00 
FOR SWEEP PATTERN 3/1 ACCEPT Ho 
Venus OF T STATISTIC 32700 


FOR SWEEP PATTERN 4/1 REJECT Hy 


VALUE OF Tl STATISTIC 34.50 


FOR SWEEP PATTERN 5/1 ACCEPT Ho 


It should be noted again that the above results were obtained 
with arbitrary sectors and arbitrary sectors values of 
detection and that the encounters were particular straight 
line encounters. In addition, the assumed characteristics 

of the sonar equipment was specific and the relation between 


Memecind=si Was GeStrretive. 


2uh 





Estimates of the probability that a target which enters 
a search pattern will be detected can be made as follows: 
Assign the values V, = Vp = Vo = 1 and V, = 90 and then 
divide the average of the values which are obtained from 
the simulation by the number of targets generated in the 
Simulation. The resulting numbers are the estimates. 
Estimates for the five sweep patterns investigated are 
given in Table III. 

TABLE III 


BpoLtintes OF THE PROBABILITY OF DETECTION 


Sweep Pattern ESt. Prob. Of Detection 
ial oka BIB 
yu 0.89 
SAA 0285 
4/1 O62 
S70 0.79 
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APPENDIX A 


Determination of False Alarm Probability. 

Since Pa and Pe apply to a resolution cell and the 
resolution cell's size in time is equal to the pulse length, 
in one hour n = 3600 Jf .1 = 36000 resolutions cells will be 
sampled. 

If independence is assumed and if X represents the 
number of false alarm in one hour period then the probabi- 


lity of less than’ two false alarms in one hour is 


af Z n k eas 
P [X<2] = L (.) Pe (1-p,) 
=O 
n n! n-1 n! 2 n-2 
Soe nem: (leek? 6 6 pee prey Pe (2 -Pp) : 
Since Pe<<l and n>>l 
n n 
P [Xs2] = (l-p,)” + n pp (1-p,) 
= (1-p,)” {l+n p,) 
Pe Pe 
= (1-np,) (1-np,;) 
= l — np,” r 
~ a 
Hence, pe = — Gime wr | Xai 9.95) that as, 
n 


Pier =. boatnen 


pe = 0.8784.10°> or 


. =5 
Pr fs 8, 


I, 
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RMAX 
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No 
RMIN 
RDA 
VSH 
VSS 
Vso 
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INCR 
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JN 
TTAB 
WEn 
MOSn 


TLXx 


APPENDIX B 


VARIABLES AND DATA IN THE PROGRAM 


Number of targets generated 

Maximum sonar search range 

Attenuation factor for propagation 
Upper limit of Sector A (starboard side) 
Upper limit of Sector A (port side) 

Seed for random number generator 
Numbers of additional targets generated 
(LIM=NSS+50) 

y distance of generation of target 
Minimum sonar detection range 

Upper range limit of Sector C 

Speed of the ship 

Speed of the submarine (target) 

DCC eons ound 

Half beam width 

Seourrnietrain Of tLhevprojector 

Number of runs 

increment in number of ping. 

Value of table for Mann-Whitney test 
Value of detection in Sector n (A,B,C,D) 
Required number of sweep in Sector n (A,B) 
Transmission loss corresponding to 0.xx prob. 


of detection 
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VARIABLES 
isyadiva 
IBEUB 
IBELO 
ALFA 

ADV 

PD 

SEL 

TL 


gull 


IGEN 


TWEn 
TTW 
NOn 


NOSn 


Angular position of projector 

Upper ray of sonar beam 

Lower ray of sonar beam 

Angular position of the target 

Advance of target in each inter ping time 
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